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We propose a field-theoretical approach to a polymer system immersed in an ideal mixture of 
clustering centers. The system contains several species of these clustering centers with different 
functionality, each of which connects a fixed number segments of the chain to each other. The 
field-theory is solved using the saddle point approximation and evaluated for dense polymer melts 
^ ■ using the Random Phase Approximation. We find a short-ranged effective inter-segment interaction 

with strength dependent on the average segment density and discuss the structure factor within this 
^\ ' approximation. We also determine the fractions of linkers of the different functionalities. 
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I. INTRODUCTION 



In this work wc consider a system where a single polymer chain is immersed in an ideal binary mixture of clustering 
centers. The study of the resulting associating polymers has a long history (see, e.g., in 0,01 and references in them). 
Models for the reversible gelation of polymers range from the consideration of pairwise associations of sticky chain 
■ segments Q to the formation of arbitrary size clusters due to association of dipolar elements in polymer chains Q ■ 
Whereas a large number of theoretical treatments model association with some form of crosslinking, i.e. linking of two 
segments only, or arbitrary size clusters, we present a treatment for a small number of species of clusters consisting 
of a fixed number of polymer segments. This falls into the so-called closed multimerization scenario [f|. 

The reversible or permanent linking of polymer chains or sections of polymers has been a topic of extensive inves- 
tigations for many decades in a wide range of systems. The general statistical physical scenario generally requires 
evaluating both the the statistical physics of the chains within a certain linked scenario as well as a summation 
over all compatible modes of linking the polymer constituents. Independent of whether the linking (crosslinking, 
aggregation, clustering, type of polycondensation, etc.) is permanent or reversible the topologically and geometrically 
i permissible combinations of linking or clustering need to be evaluated, albeit with different strategies for quenched or 
annealed situations [fj. Therefore it is natural to think in terms of the enumeration of graphs, as extensively reviewed 
by Kuchanov, et al. [31 - The ideas can then be applied to a variety of systems, such as associating telechelics 
polycondensation Q, polymers with multiply aggregating groups flO ] and general thermally reversible aggregation, 
clustering or association jll| - [l9l ] . 

There are different possibilities in which the scenario of fixed functionality clustering can be realized. One can think 
of the segments of the chain connected to particles or sidechains that assemble into structures with closed shells, akin 
to Janus particles, where it was recently shown that monodisperse ten-particle micelles, and forty-particle vesicles, 
^ ' are the thermodynamically dominant assembled structures [20h22| | . 

We are interested in the properties of a solution of such a polymer with clustering centers and in the relative 
dominance of co-existing clusters with different, yet fixed, functionalities. To this end we reformulate a ficld-thcorctical 
method originally proposed by Edwards for permanent, arbitrary- functional end- linking of chains [23| . The resultant 
field-theory is highly non-linear, but can be treated analytically and numerically, offering an additional theoretical 
tool to address questions on the formation of localized, reversible structures of polymer chains. 

As already mentioned, in order to compute the partition function or free energy one needs to evaluate the polymer 
chain conformations subject to the restrictions imposed by the functionalities of the linkers and include all possibilities 
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of linking or cluster formation. This is because we are modeling a strong type of aggregation with fixed functionalities, 
where the clusters are well-defined and local. In other words, all permissible graphs must be generated, their connec- 
tivity restrictions be imposed on the polymer chains and weighted appropriately by Boltzmann factors. Kuchanov 
et al. @, in their exposition of strategies to do this, also point out the very clear analogy between enumeration of 
spatially embedded graphs and Feynman diagrams from field theories. The obvious utility of a field-theoretical tactic 
lies in the large spectrum of available approximation techniques and graphical expansions but also in the freedom to 
choose the precise manner of implementing the additional fields. 

In this paper we introduce additional fields, with the associated functional integration, whose role is to produce the 
desirable linking, network-formation or aggregate possibilites as well as enforcing the spatial consequences of this on 
the appropriate monomers. The current approach is similar to the those used in Refs. 0, [H| and in work that can be 
seen as precursor to the current formalism 0, HE]- Whereas the specific systems investigated in Refs. [23[ and Q are 
addressed such that the ends of polymer chains or of star-like polymer units can associate, respectively, the system 
under investigation here deals with aggregation of segments of a polymer chain. We show that this system allows a 
formulation of the field theory that has the mathematically advantageous property that it is local in the introduced 
fields and these additional fields also couple to the density of monomers in a local way. 

Clearly the sum over aggregated states, by whatever method derived and approximated, would generally impose a 
complicated form onto the polymer conformational averages. This is also the case for theories with additional fields 
following integration over the fields. In approximating the functional integral over fields one expects to find nonlinear 
integral equations for these fields in the kinds of self-consistent field theory calculation that emerge in models for many 
systems (as in @, HH). However, in our formulation the saddle point equations related to the additional fields (and 
taken before integration over the polymer degrees of freedom) turn out to be only algebraic, albeit nonlinear, providing 
a significant advantage in tractability in comparison to nonlinear integral equations obtained within other strategies 
(for other systems). Analysis of the stability, etc. of the resulting theory is also relatively simpler. Completing the 
integration over the fields incorporates the clustering into the remaining weight for the conformations of the chain (also 
in a local manner) giving the "structural" contribution that is taken together with the remaining polymer-polymer 
interactions. 

Yet other path integration techniques have made use of generating functional approaches to enumerate tree-like 
configurations (26l - [28j in associating systems. Motivated by a wide range of physical scenarios under which polymer 
chains can aggregate, many different methods (mainly not field-theoretical in the sense as here) have been utilized 
in determining the contribution of certain classes of connections, ranging from summations over a subset of looped 
conformations [U, sums of tree-like graphs [29[ and trees with cycles to analyses for stickers [1, d, [3(|. Typically 
the effective polymer-polymer contributions can then be dealt with through a further self-consistent field theory (e.g. 
[3l| ) or by determining fluctuations with respect to a reference system (e.g. [lH) or through a mean-field treatment. 
In principle, before our approximations at least, the field theory introduced here is not restricted to subsets of 
connectivities or specifically cyclized conformations nor is it a priori a mean-field formalism. 

In the current calculation we have a single polymer chain that is immersed in an ideal binary mixture of pointwise 
clustering centers with different functionality (number of links that the center can have with the polymer segments). 
As a mathematical device we can think of free segments of the chain being part of clusters of functionality one, i.e. they 
cannot connect to any other segments. Moreover the system is in solution with clustering centers of functionalities a 
and b ^ a. The highly non-Gaussian field-theory resulting from the study of the model is quite complicated but can 
be approached through the saddle point approximation. Assuming the polymer to be highly dense we can then use 
the Random Phase Approximation (RPA) to describe the polymer degrees of freedom. We are then able to extract the 
local densities of segments that form part of clusters of different sizes, the effective potential based on small density 
fluctuations around a background of a given density, and the static structure factor. In the current treatment we 
develop the formalism and then investigate properties of the system in the scenario of the uniform polymer segment 
density with Gaussian fluctuations. However it is also shown where this approximation breaks down. One could 
certainly expect nonhomogenous phases to develop which can in principle also be addressed by the formalism together 
with the consideration of higher orders in density fluctuations [33l - l35j . 




We shall refer to "polymer segments" as the monomeric units of which the polymer chain is built. "Clustering 
centers" refer to point-like seeds of "clusters" of segments of the polymer chain. In other words the clustering 
centers function so as to attach to a specific number of the segments of the polymer and in so doing to localise these 
segments at a common point, binding them reversibly into a cluster. Consequently a polymer segment is also the 
basic unit to which a single attachment to a single cluster is possible. We shall deal with clustering centers of different 
functionalities which form a multi-component system with the polymer and provide a "sea" of centers to form clusters 
with the polymer segments. 

Although our field-theoretical formulation includes no precise model for the mechanism that causes clustering 
centers of a given functionality to occur, we investigate in Section l*VTl the case where the functionalities (10 and 40) 
of the clustering centers are the same as those determined for Janus particles in recent studies [22| ■ Indeed there has 
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recently been much development in the techniques for the synthesis of new patchy colloidal particles |36l - l39j . One 
particularly simple class of these anisotropic particles, called Janus particles |40l - l43l ]. seem to form mainly clusters of 
cither 10 or 40 particles. Here Monte Carlo simulations [2(| HH indicate that mainly stable micelle (10 particles) or 
vesicle (40 particles) arrangements of these particles are to be found in the vapor phase. Moreover it was found that 
the clusters behave very similarly to an ideal gas, since the particles forming the cluster tend to arrange with their 
active surfaces towards the cluster center. 

Janus chains have been suggested as potentially useful candidates for understanding interesting polymer phenomena 



44 . We will apply our formalism to the case of a dense polymer in a Janus fluid and in so doing we hope to add to 



the recent interest for Janus particles interacting with polymer chains [43j-|45j. To the best of our knowledge there 
are no results in the literature that prove the clustering in the Janus fluid in the presence of the polymer. So we will 
take as a working hypothesis the existence of such a clustering, and make the approximation of treating the Janus 
fluid as an ideal mixture of Janus partilces, micelles of Janus particles, and vesicles of Janus particles (in the spirit of 
Ref. 0). 

The paper is organized as follows: in Section |H] we describe the model we are studying and formulate the field- 
theory, in Section IIIII we perform the saddle point approximation, and discuss when we expect the approximation to 
be most accurate. Section IIVI is devoted to an investigation of dense polymer system with clustering. We use the 
Random Phase Approximation and derive an expression for the free energy density of the system and the effective 
interaction caused by the clustering centers. In Section [V] we determine the structure factor and its curvature at 
small wave vectors, in Section [VI] we finally solve the Janus case numerically, and section IVTI1 is devoted to the final 
remarks. 



II. THE MODEL 



The grand partition function for the ideal binary mixture of clustering centers of functionality a and of functionality 
b. can be written as 



e = E E E 



( Zl V) N > (z a V) N ° (z b V) N > 



iVi! N a \ N b \ 

N 1 =0N a =0N b =0 

where we also allowed for a third species of particles, of functionality one, which cannot cause aggregation. Zi are 
the usual fugacities for species i, V is the volume of the system, and Ni + N a + N b is the total number of clustering 
centers. 

We now wish to connect these to the polymer degrees of freedom and develop a partition function for the system of 
polymer together with the clustering centers. The "clusters or clustering centers" represents the free macro-particles 
making up the ideal mixture (a two component mixture with clusters of two different functionalities that are living in 
a "sea" of clusters of functionality one: the particles) in which the polymer is immersed. These macro-particles are 
made up of a fixed number of pointwise particles (i.e. they have fixed functionality) each of which is linked with one 
polymer segment. 

A suitable field-theoretic formalism was developed by Edwards [23| to describe polymer gels. Consider a field 
<j>i : M 3 — >• K then the following Wick's theorem holds (see Appendix [C]) 

I(r 1 ,r 2 ,...,r 2M ) = M j [dfa] M^i)M^) ■ ■ • M^m) e~* f dr ^ 

E S ( r h ~ T h) S ( r U -r i4 )---5(r i2M _ 1 - r hM ) , (1) 
all pairing 



where = 1,2,..., 2M and £, ^ lj for all i ^ j, and Af is the Gaussian normalization. 

If we introduce another field <j>2 , in terms of complex fields ip = (f>i + i(f>2 and Lp* = 4>\ — ifa , the following identity 
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follows (see Fig. [IJ, 

M M' 

J = Af /[#i][# 2 ] n^i( r 0+#2(r0]n^i( R i)-^2(R J )]e- /dr ^ (r) - /dr ^ (r) 

i=l j=l 
M M' 

= m' [mw] n^(^)n^( R ^ e " /drv(r) ^ (r) 

j i=i i=i 

= 5 M,Af' H 5 ( r 'i - R mi )<$(r; 2 - R„ i2 ) . ..6(ri u - R„ iJf ) , (2) 
all pairing 

where U and m, can vary over (1,2,..., M) with ^ lj and 7^ for all i ^ j . This means that each cp is 
associated with another ip* through a Dirac delta function in all possible pairwisc combinations. As shown in Fig. [TJ 
we can view the fields ip and ip* as being complementary, since the delta-function connection does not occur between 
pairs of ip or pairs of ip*. This Gaussian theory, therefore, enumerates all possible pairs of points and R_, and 
enforces this by inserting a Dirac delta. 

We consider a polymer chain consisting of TV links (see Fig. Given the Green function G(r,r ) for the segment 
of chain between two links, the conformation statistical weight of the polymer, whose conformation is described by 
the points {Ri}, is 

P({R,}) = G(Ri, R 2 )G(R 2 , R 3 ) ■ ■ ■ G(Rjv-i,Riv) • (3) 

The polymer is immersed in an ideal mixture made up of two types of cross-linked particles (see Fig. [3]) with 
functionalities a and b, respectively. We can also think of the chain segments as forming two types of clusters (closed- 
shell clusters, as for Janus particles (20[) called micelles and vesicles. What is important is the fixed functionality. 
The role of these clustering centers is to provide the links of a certain fixed functionality between polymer segments, 
thereby connecting a given number of segments of the polymer chain together. Then the partition function can be 
written as 

00 

Zn = ^2 ^Ni,N a ,N b j (4) 
N lt N a ,Nj,=0 

where TV = TVi + aTV a + 6TVb is the total number of polymer segments or of particles (since the field-theory requires 
that each segment must be paired with a particle) and 



b = Af J dRi...dR w e-' u ^"™=i i(R "- Rm) J[d>p][dp*]e-f drip{r)v * {r) x 



^(R 1 )G(R 1 ,R 2 ) V *(R 2 )G(R 2 ,R 3 ) . . . G(R N ^,IL N )p*(Il N ) x 

^ (j dr zM r)) Nl ^ (/ drz a <p°(v)Y a ^ (j Ww)"* , (5) 



and the Zi are generalized fugacities (that might also contain a multiplicity associated with the functionality). We 
have explicitly added clusters of functionality one (the single particles) here to represent the "sea" of particles in 
which the polymer and the a and b clusters are immersed. Clearly clusters of size one simply link to a single polymer 
segment and therefore do not cause association as the a and b clusters do. Note that we are not interested in describing 
the precise model for the mechanism by which the clusters of a given functionality are formed from the aggregation 
of particles (more on this in Section fVij) ; we just assume that this aggregation process takes place. At this level 
of description the physics of the precise mechanism for the multimerization resides in the fugacities. In our present 
formalism the clustering centers and the particles are pointlikc. In principle it is also possible to extend the current 
formalism to model clusters of finite extension. We have also added an excluded volume term v, to the polymer chain, 
with the dimensions of a volume. 

It is possible to use the formalism without necessarily introducing the essentially inert clusters of functionality one, 
which turns out to be the fugacity z\ = 1 case of the equations derived below. Appendix [Dl shows the details. We 
continue with the slightly more general formalism here, noting that Z\ — > 1 will show no effects due to the addition of 
these convenient clustering centers. 
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We then find in a short hand notation, and neglecting for the time being the excluded volume term, 
Z N = Af J[di P ][dip*]{l[dR}{l[G}r l N e X p(- J dv <p(r)<p* (r)+ 

J drp(r)ln(tp*(r)/r))+zi J drip(r) + z a J drip a (r)+z b J drip b (r)j (6) 
= Af J [dtp][dtp*] {lI dR } {[[G}eM^,^}} , (7) 



where we introduced the microscopic density of polymer links p(r) = S(r — R^) and rj is an arbitrary constant 

with the dimensions of a length to the power —3/2. In the rest of the paper we will measure lengths in units of r/~ 2 ^. 
A natural choice would be r\ = £~ 3 / 2 : with I the Kuhn length of the polymer segment. 



A. A simple example 

To clarify our formalism we here consider the simple example of a polymer, with 4 polymer segments, interacting 
with two clustering centers of functionality a — 2. Using the properties of the Gaussian chains, the partition function, 
neglecting the excluded volume term, is written as 



z = u j [mw\ {n dRi | {n**} |n G ( R i' R i+i)| 

jn^(R0jz 2 V(riV 2 (r 2 ) 
= Mz\ J dridr 2 G(r fl , r i2 )G(r i2 , n 3 )G(r i3 , r i4 ) . 



- / dr ip(r)ip* (r) 



(8) 



where we used Eq. ([2]) and the prime on the last summation symbol indicates that we have to sum over all possible 
ways of assigning to the indexes (ii, 12, 13, i&) the set (1, 1, 2, 2). We then sec that the result is given by 



Nz 



4 J dridr 2 G(ri,ri)G(ri,r2)G(r 2 ,r 2 ) +4 / dridr 2 G(n 3 r 2 )G(r 2 , r 2 )G(r 2 , n) + 



4 J dridr 2 G(ri,r 2 )G(r 2 ,ri)G(ri,r 2 ) 



— °"oO 




(9) 



where lines represent polymer segments and black circles represent integration. 
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III. THE SADDLE POINT APPROXIMATION 



We look at the saddle point equations for the fields tp and p* for any arbitrary yet fixed conformation of the polymer 
p = p(r) > Vr G K 3 . The saddle point approximation becomes necessary because the field-theoretical formulation 
of the system, Eq. ([5]), is certainly highly non-Gaussian. We proceed to calculate the saddle point and note that, for 
fixed polymer conformation, the resulting condition is a set of algebraic equations, i.e. the saddle point depends only 
on the local density. We shall also see that there exists always exactly one saddle point, corresponding to the physical 
solution of the equations. 

We need to determine the critical point (ip, <p*), which amounts to solving 



ST 



5<p(r) 



<p,<p* 



= 



ST 



5(p*(r) 



= 



if, if* 



(10) 



These can be combined to give 



- <p* + zi + az a p a 1 + bz b Cp b 1 = 
-tp + p/tp* = 

p/tp + zi + azap^ 1 + bz h p b ^ x = . 



(11) 
(12) 

(13) 



A. Properties of the solutions of the saddle point approximation 



From Eq. (|13p we can write 

-p+m = o, (14) 

where f(x) = Z\x + az a x a + bz b x b . We can then show that there exists at least one positive solution of Eq. (fTS")) for 
p > 0. Indeed when <p — > then — p + f((p) = —p < 0, whereas when ip > pj z\ then — p + f(<p) > 0. Consequently 
— p + f((p) goes through a zero when ip G [0, oo). And it is easy to show that there is exactly one positive solution 
(p > for any p > 0. f'(0) = Z\ + a^Zap^ 1 + b 2 zi,p b ^ 1 so f'(0) > for all p > 0. Since the derivative is strictly 
positive and the polynomial goes from negative to positive, there exists exactly one positive root (p for any p > 0. As 
we show in Section HVBI the local cluster densities of the saddle point are proportional to z\(p, az a (p a , and bzbp b , such 
that Cp must be non-negative. 



B. Solutions for ip and <p* 



Consider the particular case a = 2, b — 0, corresponding to reversible crosslinking, then from Eq. fjl3|) we find for 
the critical point 



_ ± = z z 1 ±^BEe , (15) 

— Z\ ± \J Z{ + &Z2P 



which rules out the p- solution, as the logarithm of p*_ is not well defined. 



C. Expansion of T around the critical point 



We can expand the function T around the critical point to second order in the fluctuations Atp = tp— <p = A</>i+iA02 
and Aip* = tp* — p* = A^i — iAfo, 

T[<p + A(p,(p* + Atp*] = T[p,p*} + 

Ifdnj dr 2 ( A^(rx) A0 2 (n) ) . F 2 (rx,r 2 ) • ( 
third order terms , (17) 
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where -F-^ri, 1*2) is the following 2x2 matrix 

F 2 (n,r 2 ) = 



& 2 r 

<50i(ri)50i(r2) 
8 2 F 



(502(ri)(50i(r2) 



S 2 JF 
S(t> 1 (r 1 )S<t>2{r 2 ) 

S 2 T 



/>2(ri)502(r2) 



<P,<P 



We have 



(5^i(ri)^i(r 2 ] 



-2 - + a(a - \)z a (p a - 2 + b(b - \)z b (p 



(50i(ri)^ 2 (r 2 ] 
5 2 F 



— + a(a- l)z a Cp a - 2 + b(b - l)z b p 



b-2 



S<j} 2 (ri)S(j)2 (r 2 ) 



-2 + ^ - a(a - l)z a $ a - 2 - b(b - l)z b (p 



b - 2 6(n-r 2 ) , 
5(n - r 2 ) , 
! <y(ri-r a ) , 



We can then introduce the matrix M. [p] as 



-2 --4. 



+ a(a - \)z a Cp a - 2 + 6(6 - l)z b $ b ~ 2 i + a(a - X)z a (p a ~ 2 + 6(6 - l)^ 6 " 2 ] 



+ 0(0 - l)^' 2 + 6(6 - l)z b £ b 



-2 + -4j - a(a - \)z a (f> a ~ 2 - 6(6 - l)z b (p b 



Then the partition function can be rewritten as 



where 



and 



Z N ^N' J {J] dR] {[] G} exp{M[p]} , 
M[p] = T\(p, dr In {det M[p}} , 



det M[p] = 4[l + z a a(a-l)(p a /p + z b b(b-l)<p b /p] 



(18) 



(19) 
(20) 
(21) 

(22) 

(23) 
(24) 
(25) 



For finite p the determinant is always positive with positive real parts of the eigenvalues, if tp > 0. Hence the saddle 
point is stable. 

As the density p increases we note that the free energy calculation due to the saddle point F[ip, (p*\ grows at least 
as p, but that the fluctuation contribution — \ J dr In {det strives to a constant. Therefore at sufficiently high 

density we expect the relative contribution of the fluctuations to become negligible. 

In the particular case a = 2, b = we find 



detM±[p] = 4(1 + 2*2^/^) . 



(26) 



IV. THE RANDOM PHASE APPROXIMATION 



In this section we express the polymer in terms of a segment density via the so-called random phase approximation. 
We restrict ourselves to the situations where the polymer segments are presumed to be distributed homogeneously with 
Gaussian density fluctuations around this value. For sufficiently dense systems this has been treated as a reasonable 
approximation [6j in permanently networked systems. We do not investigate the cases of possible inhomogeneous 
phases. In principle we would then need to expand our results in the preceding sections to hi gher orders in the 
density fluctuations [33T - [35| or attempt to express our results in terms of more complex quantities |46l. |47| . However, 
we do investgate where our RPA with the homogeneity assumption fails as one indicator for possibly different phase 
behaviour in the system. 
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A. Basic formulation 



Clearly our clustering formalism produces a significantly nontrivial density dependence. Presuming that our system 
behaves like a highly dense polymer melt, where the fluctuations of p = p + Ap are small, one can use the following 
Random Phase Approximation [35j (RPA) (see Appendix [Bj , 

/ {n rfR }{ll G }--- = Af" J[dAp}e X p(-l JdrJ dv Ap{v) S « 1{]V v r l) A P (r ) 

= N" J [dAp] exp ( -II £ Ap(k)^i^Ap(-k) j . . . , (27) 



where we denoted with a tilde the Fourier transform [48j and with a hat an inverse Fourier transform. We note that 
this type of approach is not atypical in calculations for quenched gels |4{| . 

Expanding to second order in the density fluctuations we find, from Eq. (|2~i)) and Eqs. pT|) - ([T2l . 

M[p + Ap] = AV + B J drAp(r) + C J dr[Ap(r)} 2 + ... (28) 
= AV + C J dr[Ap(r)} 2 + ... , 

where we used the fact that J dr Ap(r) = 0. V is the volume of the box, A, B, and C are given functions of Zi, 22, and 
p and for the particular case a = 2, b = can be found in Appendix [A] Notice that in this case A_ is not defined for 
any values of the average density so only the A + , B + , and C+ solution is physically meaningful, i.e. they correspond 
to the expected positive (p solution (see Section IIII Ap . 

We then obtain the following approximation for the partition function of Eq. (|23[) 

Z w «AA"'e AV '-5/*in(J?/v-2C) ; (29) 

where i^" 1 is the operator whose r,r' component is given by <S ( 7 1 (l r — r 'D- 

In terms of the free energy density /3f = — ln(Z^)/V we find, in the thermodynamic limit (V — > 00 with p = N/V 
constant), 

/?/ = -A-\ H-2C) . (30) 



B. Local clustered segment densities 

Following the usual method for grand-canonical ensemble it is possible to calculate the local densities of segments 
that form part of different sizes of clusters. We then compute 

zid\nZ N 

ni = iv^T' (31) 

for segments of the chain that are not cross-linked, and 



_ z x 8\nZ N 

n dz x ' (6Z) 



for the density of segments part of clusters of size x. 

Neglecting the logarithmic corrections due to the quadratic fluctuations, we have 



hiZ N = VA = M[p] . (33) 



Within the saddle point approximation we split the contributions into parts due to the saddle point above and due 
to the quadratic fluctuations (notice that this analysis holds also locally at the level of the partition function for an 
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arbitrary p = p(r) profile of polymer chain density) 

z x dT[(p, <f*]\p=p _z x V dlndet M[p] 



~N dz x X ~2N dz x 

= nf + n« , (34) 

where x = 1, x = a, or x = b. We then find 

n s x = ^f. (35) 

Since n x have to be real and non-negative, in the saddle point approximation, also the solutions (p have to be real 
and non-negative. Immediately from the saddle point equation (|13[) follows that 

nf + n s a + nf = 1 , (36) 

must hold generally. We also find after some algebra 

nf + n® + = . (37) 

This means that the saddle point approximation conserves the total number of segments for any density. Consequently, 
any average over these density dependent expressions, irrespective of the approximation, must satisfy the conservation. 
However we note that the conservation laws (|3"6"|) and ([3T)l do not prevent possibly negative nf + n® which can arise 
in the region where the fluctuation part is not sufficiently smaller than the saddle point. As the validity of the saddle 
point improves with density, also this possibility disappears. 
For the special case a = 2, b = we have 

2Zl (38) 



z x + \Jz\ + 8z 2 p ' 



2p{z\ + 8z 2 p) 



as a consequence we see that the fraction of monomers not in a cross-link decreases with the density. 

C. The effective potential 

Upon integrating over the degrees of freedom associated with the clustering centers the remaining integral in the 
partition function is that over the polymer density degrees of freedom (in the RPA). This permits us to interpret 
the effective interaction between polymer segments as caused by the clustering. It consists of the typically attractive 
contribution to the polymer-polymer quadratic density fluctuations from the aggregating fields and any direct polymer- 
polymer interaction (such as excluded volume interactions). 

From Eq. (|23p and (|2"T)) - (|28[) we can rewrite the partition function as 

Z N « Af"' f [dAp}e~ ^ dr f dr ' Ap(r) S ° 1 'v r ' r - A ^ r ' ) e AV e -h$ d *S dr ' Ap(r) W(r-r' ) Ap(r ) ^ 



where the effective potential between the polymer segments, W, is given by 



W{r - r') 



5 2 M[p] 



= -2C( Zl ,z a ,z b ;p)S(v-r') . (41) 

P(r)=P 



Sp(r)Sp(r') 

We can then split the contribution from the saddle point and the quadratic contribution of Eq. (|24|) and write 



C s 
Co 



C = C s + Cq , (42) 
1 d 2 f s &,p) 



2 dp 2 

1 d 2 fQ(^ P ) 

2 dp 2 



(43) 
(44) 
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where 



I S = ~P + pki(p/tp) + znp + z a ip a + Zbip b , 
f Q = ~\ ln[4(l + o(o -l)z a p a /p + 6(6 - l)<p b /p)] • 

Now, using Eq. (fi"3]l. we find df s /dp = hip — In 9? and , using the property d(p/dp = l/(zi 



a 2 z a ip a 



1 + b 2 z b <p b - 



(45) 
(46) 



follows d 2 f /dp 2 = 1/p — l/(z\(p + z a a 2 (p a + Zbb 2 (p b ). Let us assume for definiteness that b > a. Then when p is very 
small z\(p « p and d 2 f s /dp 2 « a(a — l)(z a /z")p a_2 /[l — a(z a /z°)p a_1 ], while when p is very large Zbb(p b w p so that 
d 2 f s /dp 2 « (6 — I) /(bp). Moreover we find in the large p limit that d 2 f®/dp 2 behaves at least as l/p 2 . 

We remark that in the small p limit in the a = 2, b = case d 2 f® /dp 2 ~ I0z 2 /zf whereas in the a = 10, 6 = 40 
case d 2 f Q /dp 2 w -3240z a p 7 /z 1 10 . However we guard against interpreting this as a repulsive interaction as the saddle 
point approximation to our field-theory is not expected to be accurate at small densities. This repulsive contribution 
in the small density limit for the effective potential of the Janus case (see Section IVI[) explains the fact that here the 
RPA can be valid (see Section |IVD[) even if we do not add any excluded volume interaction to the polymer. 
For the simple case a = 2, b = we then find that C+ = C5 + Cq where 



C s 




Here we also find in the small p limit 



Cc 



Z'2 



fe 2 



l±-^lp + 0(p 2 ) 



5z 2 



0(P 2 ) 



(47) 
(48) 

(49) 
(50) 



which tells us that the energy to cross-link goes to a constant proportional to Z2- 
In the large p limit 



Cs 
Cq 



1 Zl 
4p 8V2i^P 3/2 
3zi 



128V2z|p 5 / 2 

+ o(p- 3 ) , 



+ 0(p" 3 ) 



(51) 
(52) 



which tells us that the energy to cross- link goes to zero as l/p, in accord with the fact that we are in a dense system. 

The effective potential calculated here, based on small density fluctuations around a background of a given density, 
is dependent on the average density. As expected the clustering produces a local attractive interaction. 

However, it is interesting to note that the strength of this interaction decreases with average density. We attribute 
this to the fact that the fraction of free segments (i.e. those in clusters of size 1) decreases with the average density 
according to Eq. (|38[) . Therefore, for large p, the number of additional free segments gained by increasing the density 
from p to p + A is proportional to p~ x / 2 leading to a pairwise contribution A 2 /p as found in Eq. (|5T|) . 



D. Validity of RPA 



The RPA is based on a homogeneity assumption which no longer holds when the RPA itself predicts overly large 
fluctuations. In order to obtain Eq. (|2"9"|) we must have that S () ' 1 (k)/V — 2C is a strictly positive function for all values 
of the wave vector k. Since SQ 1 (k) is a monotonically increasing function of k, the RPA will be valid as long as 



C< 



Sp HQ) 

2V 



2Vp 2 



(53) 



In the thermodynamic limit one would require that C < 0, for the validity of RPA. 
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We can then extend the region of the validity of RPA by adding an excluded volume effect [50| to the polymer 
which amounts to taking M[p] — > M[p] — v J dr p 2 (r) with v a positive constant with the dimensions of a volume. We 
will then have 

A-tfv , B -> B -2pv , C -> C - u , (54) 

and the validity of RPA becomes C < v. 

For the a = 2, b = case we have that C+ is always positive so the RPA cannot be applied without the excluded 

volume interaction. As a matter of fact we have lhxip^oo C+ = and limp^o C+ = z i{z\ + 5z2)/zf and C+ is a 

monotonically decreasing function of p. So by choosing v any arbitrarily small positive constant we are able to extend 

the range of validity of RPA to arbitrarily large densities. 

In this case choosing z\ = 1,22 = e 2 ^, and v = 1 the validity domain in the phase diagram is determined in Fig. 2] 
In Fig. [5]we show the behavior of the free energy density as a function of density in the case a — 2, b = 0. Here we 

choose Z\ = 1, Z2 = e 2 @ were (3 = with Boltzmann constant and T is the temperature and v = 1. 



V. THE STATIC STRUCTURE FACTOR 

The Fourier transform of the pair correlation function is defined as (5l| 

g(k) = I(p(k)p(-k)> = i<A/S(k)Ap(-k)> , MO. (55) 

The quantity g(k) can be measured experimentally by light scattering. Moreover one can extract some important 
information on the polymer properties from the small k = |k| behavior: 

g(k) = .g(0) (l - y R 2 g + . . .) , (56) 
where R g is the radius of gyration of the polymer, namely 

N N 

^ = ^ E E « R » - R ™) 2 > • (57) 

71— 1 m. — 1 

Now using the result from the RPA we find 

N~g{k) = V[S^ l {k)/V ~2C]~\ (58) 

which when C — — v, agrees with Edwards' result (HoL l52| for polymer chains with only excluded volume interactions 
[S3j . Notice that the effective potential C is in general a function of zi, z a , and z&. As shown in Section HV CI C + v 
tends to be positive (attractive interaction between polymer segments) in the presence of clustering centers. So we 
expect there to be a regime of density for which there is a balance between the repulsion due to the excluded volume 
effect and the attraction due to clustering making C « 0. In such case our result reproduces the one for the ideal 
chain (see section 1.2.3 in Ref. (5l|). 
In the small k limit we find 



Ng(k) = 



l/p 2 -2CV 



k 2 
3 ^ 



(59) 



where p = N/V is the average polymer segment density for a single long polymer chain and the "curvature" of the 
structure factor at k = is 



2 SpV 

? = 6 - Up^CV ' (60) 



where i is the Kuhn length of the polymer. And in the thermodynamic limit 

(«/<)■-» -5^5- ( 61 ) 

So that at large polymer densities the curvature tends to a constant (we note that C now includes the excluded volume 
as in Eq. ([M]) ). 
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We also find, in the thermodynamic limit, the following expression for the structure factor 

^ - W^mc ■ (62) 

Notice that in the absence of the effective interaction (C = 0) the structure factor diverges at k = 0. 

In the a — 2, b — case, at constant V, in the small N limit we find the free polymer result £ 2 = N£ 2 /6 + 0(N 3 ). 
In the large N limit we find 

Given the densest possible filling N/V ~ 1/v the curvature tends to a constant. 

VI. THE JANUS CASE 

Although our field-theoretical formulation includes no precise model for the mechanism that causes clustering 
centers of a given functionality to occur, we investigate here the case where the functionalities (10 and 40) of the 
clustering centers are the same as those determined for Janus particles in recent studies. Indeed there has recently 
been much development in the techniques for the synthesis of new patchy colloidal particles [36l - l39| . One particularly 
simple class of these anisotropic particles, called Janus particles |40l443| . seem to form mainly clusters of either 10 or 
40 particles. Here Monte Carlo simulations [2(3, Hl[ indicate that mainly stable micelle (10 particles) or vesicle (40 
particles) arrangements of these particles are to be found in the vapor phase. Moreover it was found that the clusters 
behave very similarly to an ideal gas, since the particles forming the cluster tend to arrange with their active surfaces 
towards the cluster center. 

Janus chains have been suggested as potentially useful candidates for understanding interesting polymer phenomena 
[44j . We will apply our formalism to the case of a dense polymer in a Janus fluid and in so doing we hope to add to 
the recent interest for Janus particles interacting with polymer chains (43l - |45j . To the best of our knowledge there is 
no results in the literature that proves the clustering in the Janus fluid in the presence of the polymer. So we will take 
as a working hypothesis the existence of such a clustering. And make the approximation of treating the backbone 
units of the polymer (the Janus particles) as an ideal fluid. 

Given the general setting described above we can apply our theoretical model to a polymer in a Janus fluid. By 
this we think of chain segments only clustering to form limited closed shell conformations, i.e. micelles and vesicles. 
As mentioned before, we however do not consider the nature of spacial extent of the clustering in detail and simply 
presume that it still occurs in the same way as if the Janus particles where not connected to the polymer. 

In the Janus case we have to choose a = 10 and b = 40 (see Fig. [3]). We then find for the determination of the 
critical point an algebraic equation of degree 40, Eq. (|13[) . As expected, this equation has just one solution for which 
A (from Eq. (|28p) is real and non- negative. 

We can see the generalized fugacities defined as Z{ oc exp(— ftm + fiHi) for i = 1,10,40, where m is the average 
internal energies of the cluster of i Janus particles and /i,; is the chemical potential of this cluster species. It is moreover 
reasonable to take /jj w /j independent of i (fj, being the chemical potential of the vapor phase of the Janus fluid) so 
that we get z% oc exp(— /3ui). 

Choosing z\ = l,2io = e 10 ' 3 , and 2:40 = e 40 * 3 we find that at small densities (where the theory is expected to be 
not good) the effective potential is repulsive (due to the quadratic fluctuations in the theory) even without adding an 
excluded volume to the polymer (see Section IIV CP . The range of validity of RPA in the phase diagram is shown in 
Fig. [5] In the same Figure we show the behavior of the free energy density which clearly shows the signature for the 
breakdown of the RPA theory at high density. 

Notice that, also in this case, lim /5 _j. 00 C = so that by adding a small excluded volume will allow to reach the 
high densities domain with RPA. At a fixed temperature C, as a function of density, has a global maximum, so that 
choosing the excluded volume v bigger than this value, the RPA can be made valid at any density (see Fig. [7]). 
Moreover we expect the theory to give consistent results in the high density regime. 

Our choice for the fugacities is justified a posteriori since for /3 < 1 we are in the high temperature regime of the 
Janus vapor poj where the internal energy of a cluster of i Janus particles is with a good approximation given by 
— (i — 1) rj — i (corresponding to a completely stretched cluster). 

Since Zn is a grand-canonical in the clusters of Janus particles, we can take derivatives with respect to the 
generalized fugacity Zi to determine the concentration rz.j of clusters of i Janus particles as follows 
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where lnZ N = V[A + ln(-2C)/2]. 

A graph of the concentrations as a function of the average polymer density is shown in Fig. [5J From Fig. [5] one can 
see the difference between n x — nf and n x = nf + for x = 1, 10, and 40 and 0=1. Note that the conservation 
of particles n\ + nio + = 1 is exactly satisfied at all densities and temperatures in both cases. In Fig. Owe show 
the dependence of the concentrations from the temperature. We thus would say that at sufficiently high densities the 
vesicles appear and as a consequence the micelles are reduced. 

We conclude that this suggests strong dominance of non-clustering at low densities. As the density is increased 
smaller clusters and eventually larger clusters dominate the linking behavior. 



VII. CONCLUSIONS 



In this work we have studied and developed a field-theoretical formalism for a polymer immersed in an ideal mixture 
of clustering centers. These centers cause clustering of either a particles or b particles, i.e. clusters of either species 
are monodisperse. The field-theory couples fields associated with stickers to the polymer chain density and provides 
a formally exact expression for the partition function (canonical in polymer and grand-canonical for the clustering 
centers). We showed that it is possible to compute quantities using the nonlinear theory by means of a saddle point 
approximation and we argue that the approximation improves as the density of the polymer chain increases. The 
current system and the choice of implementation of additional fields enabled us to derive saddle-point equations that 
are simpler than those that arise in some other formalisms by not requiring the solution of nonlinear integral equations. 
The benefit of the local saddle-point equations is that they also enable a relatively simple analysis of the stability and 
applicability considerations of the theory. 

For a homogeneous, dense polymer system, we computed the effective interaction potential (up to quadratic density 
fluctuations) and computed properties of the structure factor within the random phase approximation. As expected, 
the addition of an excluded volume interaction will compensate for the attraction due to aggregation effects and 
extend the validity of the RPA. The effective pairwise potential obtained in this approximation has interesting, 
nontrivial density dependence. Another clear consequence of increasing chain density is the growing dominance of the 
higher-functional clustering centers. 

The nature of the clustering process is definitely of importance in aggregating polymer systems. (Recently, a theory 
for cluster formation in homopolymer melts was introduced by Semenov fl7j|.) One motivation for our study is the 
closed multimerization scenario suggested by particles in a Janus fluid, where micelles (a = 10) and vesicles (b = 40) 
are known to occur [20T l2l| . Future work will focus the attention on the stability of such Janus- type multimers when 
connected to a polymer with more detailed models of the cluster itself. 



Appendix A: Coefficients A, B, and C for the a — 2, b = case 



The A, B, and C coefficients arc given by 



A± = -f— { -z{ ± zi \lz\ + 8pz 2 + 8pz 2 \n 

8z 2 



Apz 2 



-zi ± ^Jz\ + 8pz 2 



4z 2 



p + In 1 ± 



Zl 



\Jz\ + 8pz 2 



lni 



B± 



2z 2 



zf + 8pz 2 4p 



-h it 



zi 



( z ( + 8pz 2 y 



+ 



-1 ± 



+ 8pz 2 1 

Zl 



In 



4pz 2 



y/zf + 8pz 2 



P 2 T 



zi ± \Jz\+ 8pz 2 J 
2z x {z\ - 2z 2 + 8pz 2 ) 



(z2 + 8pz 2 ) 3 / 2 



(Al) 
(A2) 
(A3) 
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Appendix B: The Random Phase Approximation 

For the polymer chain, with a Kuhn length £, we can write 

/{n^}{iM= 

J [dp] J [dC] J [dR] e -*/o^^ 2 ( S ) e ^*C(r)[p(r)-/ ^ 5 ( r -R( S ))] = 



[dp] [d(] /[dRje-AJ'o dsR2 ( s ) e l / dr C(r)p(r 



• / yC(r)«(r-R(*))- 



1 



ds 



- dr _ C (r)5(r-R( S )) dr T C(r)i(r - R(* )) 



(Bl) 



where L = N£ is the total polymer length and the dot denotes a derivative with respect to s. 

Now the first term gives just a normalization constant TV. To calculate the second term we introduce the polymer 
center of mass R so that R(s) = Ro + AR(s) and write 

d(] /[dR]e-*J"o i ^H- 2 W e */*CW/>W f dr / L ^C(r)<5(r-R(.s)) = 



idC] / ^R] e"*^ ^R 2 M e4 /drC(r)p(r) / «fr y CW^ ^ e ^r- lk R( S ) 



J KK /drC(r)p(r) A^ /" drC(r) , 



ikr e -ik(R +AR(s)) 



(B2) 



where V is the volume of the box and (5k. o is the Kronecker delta. 
The third term gives 



|[ dCK /drC( r )p(r) y [dAR] J dRoe -4/,f (i .AR» W _^ J-J- | dr jf L & | ^ ^ ^ 



k k' 



C(r)C(r )e Ikr e 



'\„ikr ik r i(k+k )R ikAR(s) ik AR(s ) 



-R(s )) 



where we denoted with the average 



■>o 



/[dR]e-4/o dsR2 «[...] 
/[dRJe-A/o^^ 2 ^) 



(B3) 



(B4) 



and with the tilde the Fourier transform [481 ] . 

Now the average ^ e 4k ( R ( s )- R ( s ))\ can be easily calculated by discrctizing the polymer and integrating over the 
bond vectors = Rj+i — R», as follows 



,»k(R(a)-R(s )) 



/ IL db t e-^ E - b 2 elk ( bl +b 2 +...+b„ 



fUidb.. 



e 2<- 



= ( e e 



where n = \s — s \/£. Since 



ds / ds' e - als - s 1 = 
Jo 



I 2(aL-l + e- ai 



(B5) 



(B6) 
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we can introduce the function 



So(k) 



72(k 2 £L/6-l + e-^) 
V 2 kH 4 



(B7) 



with 5 (0) = (N/V) 2 = p 2 . 

Then the expression we started with in Eq. (|B1[) can be rewritten, omitting the functional integral over the density 
collective variable, as 



TV / [d(]e^ 



£ k C(k)p(-k) 



1 4 C ~ (0) - W E C(k)yS (fc)C(-k) 



We can now reconstruct the exponential to obtain, 

"[^] e ^^f( k )K- k ) e -wE.cC(k)^ o (fc)C(-k)-i^c(0) 



W' e - W E k ln(S (*W e -W E k Ap(k)^L_i^ A p(-k) _ 



Here Ap(k) = /5(k) - pF5 kj0 - 



(B8) 



(B9) 



Appendix C: The Gaussian distribution 

The Gaussian distribution function for a set of real variables x\, X2, ■ ■ ■ , Xjv is defined as 

*(a;i,x 2 , . . -,x N ) = Cexp -- ^ A nm x n x, 



(CI) 



where A rmi is a symmetric positive definite matrix and C is a normalization constant given by the requirement 

J — 00 J — oc J- in « 

Let (. . .) be the average of the distribution function of Eq. fjCll) , 



/OO /"OO 
• • • / ~[dx 
-oc J —00 



then, it can be proved [50|, that 

In general we have the following formula (Wick's theorem) 

(x ni X n2 . . . X n2p ) ^ ^ (^mi ) (^13 ) 

all pairing 



(C2) 



(C3) 



(C4) 



If the subscript n ofx n is regarded as a continuous variable, the set of points (x±, X2, ■ • ■ , xn) represents a continuous 
function, and the integral over the set (xi, X2, ■ ■ ■ ,xn) reduces to the integration over all the function, and it is called 
the functional integral. It is denoted by the symbol [da;], i.e. J JJ n dx n ...—>■ J[dx) 

Consider now the following Gaussian distribution functional 



*[</>] = Cexp 



dx<p 2 (x) 



where <j> is a real function, then using the continuous limit of Eq. (|C3|) we find 

(4>{x)4>{x')) = 6(x-x') . 

where S is the Dirac delta function. 



(C5) 



(C6) 
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If now <p — 0i + i<f>2 is a complex function we consider the Gaussian distribution functional 



Now we find from Eq. (IC6I) 



dx4>(x)4>* (x) 



(C7) 



(0(x)0V)) - (0i(x)^i(x') + ( / ) 2(x)0 2 (x')) 



and 



i<J(a: - a') + ^(a? - a') = <S(a; - a/) , (C8) 



(<f>(x)(j>(x')) = (Mx)Mx')-Mx)fo(x') + i<l>i(x)Mx') + iMx)Mx')) = 

±6( x -x')-±6(x-x')=0. (C9) 



Appendix D: Field theory without 1-clusters 

An alternative way to formulate the clustering, without the use of clusters of size 1 is presented below. We shall 
show that a simple mapping reduces again to a special case of Eq. 

Consider a system in which we have only clusters of sizes a and b but no "inert" clusters of size 1. As explained in 
Section [XT] the functional integration over the fields tp and tp* requires matching each ip of a clustering center with a 
ip* on the polymer, permitting no unmatched tp and tp* pairs. Since size a and b clusters do not necessarily attach to 
each potential site on the polymer, all possible attachment sites have to be generated. The product 

JV 



1(1 + 



i=i 



produces all equally weighted possibilities of the attching to the the sites {-R;} , Vi £ {1, . . . , N} of a given polymer 
configuration. 

The analog to Eq. ([5} then becomes 

z ' Na ,N b = W dRi...dR w e-" E "™= l5(R "- Rm) f [dip][d<p*]e-^ dTV ^^ T) x (Dl) 

(1 + ^(Ri)) G(Ri,R 2 ) (1 + v>*(Ra)) G(R 2 , R 3 ) . . . G(K S -i, Riv) (1 + V*(Rjv)) x 
Iff \ N « i / /■ \ Ar, > 



; drz a ^(r)j —l^J d vz bV b (r)j , (D2) 
leading by the same procedure as described in Section [TT| to analog of Eq. ([5]) 

Z' N = J\f J[dcp][dp*]{l[<m}{l[G}d N exp(- J drcp(v)<p*(r)+ 

drp(r)\n(l + tp*(r)/d)+z a [ drtp a (r)+z b [ drtp b (r)) . (D3) 



We see that the trivial transformation tp* — >• <p* — 1 in Eq. (|D3|) above leads to the orgininal field-theoretic equation 
in the main text, Eq. (j6|), with z\ 1. For this reason we treat the marginally more general case in this paper. 
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FIG. 1: A schematic representation of the role of the field-theory. The fields ip and tp* are depicted as functions of spatial 
variables. Multiplication by exp(— f ipip*) and subsequent functional integration enforces the linking of the spatial coordinates 
between pairs of ip and tp* (in all possible ways). 
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FIG. 2: Shows the polymer made up of N equispaced links that are susceptible to being linked into clusters. 




FIG. 3: Shows the clusters of Janus particles: the micelles are made of a = 10 links whereas the vesicles of b = 40 links. 
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FIG. 4: (Color online) Shows the RPA validity region of the phase diagram, in the a = 2, b = case, for Z\ =\,z% = e 2lS and 
v = 1. 
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FIG. 5: (Color online) Shows the free energy density as a function of the average density in the a — 2, b — case, for 
zi = 1, Z2 = e 2li and v = 1. 
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FIG. 6: (Color online) The upper panel shows the RPA validity region of the phase diagram, in the Janus case, for zi = 1, zio = 
e 10 ' 3 , Z40 = e 40 ^, in the absence of any excluded volume effect. At j3 = the C = equation has solution p ~ 0.647933 
Note that the validity region is in the small density region, where the contribution from the quadratic fluctuations of the theory 
dominates, and the whole theory is expected to be less significant. The lower panel shows the free energy density as a function 
of the average density. The rapid increase at high density is indicative of the limit of the RPA applicability. 
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FIG. 8: (Color online) In the Janus case for z\ = l,zio = e 10,3 ,Z4o = e 408 , v = 15, and /3 = 1 shows the concentrations of 
clusters of 1, 10, and 40 particles as a function of the density. 
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and v = 15, shows the concentrations of clusters of 



FIG. 9: (Color online) In the Janus case for z\ = 1, zio = e 1013 , Z40 

1 (nf ), 10 (nfo), and 40 (nf ) particles as a function of the density when we do not use the logarithmic correction in Eq. (|24|l 
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Figures captions 

Figure 1: A schematic representation of the role of the field-theory. The fields ip and p* are depicted as 
functions of spatial variables. Multiplication by exp(— J (pip*) and subsequent functional integration enforces the 
linking of the spatial coordinates between pairs of ip and ip* (in all possible ways) . 

Figure 2: Shows the polymer made up of N equispaced links. 

Figure 3: Shows the clusters of Janus particles: the micelles are made of a = 10 links whereas the vesicles 
of b = 40 links. 

Figure 4: (Color online) Shows the RPA validity region of the phase diagram, in the a = 2, b = case, for 
Z\ = 1, 22 = e 2/3 and v = 1. 

Figure 5: (Color online) Shows the free energy density as a function of the average density in the a — 2, b = case, 
for z\ = 1, z% = e 213 and v = 1. 

Figure 6: (Color online) The upper panel shows the RPA validity region of the phase diagram, in the Janus 
case, for z\ = 1, z\q = e 10/3 , Z40 = e 40/3 , in the absence of any excluded volume effect. At j3 = the C = equation 
has solution p ~ 0.647933 . . .. Note that the validity region is in the small density region, where the contribution 
from the quadratic fluctuations of the theory dominates, and the whole theory is expected to be less significant. The 
lower panel shows the free energy density as a function of the average density. The rapid increase at high density is 
indicative of the limit of the RPA applicability. 

Figure 7: (Color online) Shows the free energy density as a function of the average density in the Janus 
case, for z\ = 1, Z\q = e 10 * 3 , z 40 = e 40/3 , and v = 15. 

Figure 8: (Color online) In the Janus case for z\ = 1,-zio = e 10 @,Z4o = e 40 ' 3 , v = 15, and (3 = 1 shows the 
concentrations of clusters of 1, 10, and 40 particles as a function of the density. 

Figure 9: (Color online) In the Janus case for z% = l,z%o = e 1013 ^^ = e 40 ^, and v = 15, shows the concentrations 
of clusters of 1 (nf), 10 (nf ), and 40 (nf ) particles as a function of the density when we do not use the logarithmic 
correction in Eq. (j2~4")l . 



